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We investigate the conditions leading to polaron formation for a single electron interacting with 
dispersionless optical phonons within the Holstein and the Su-Schrieffer-Heeger models. Both from 
analytic perturbation theory and exact numerical diagonalization of small clusters, we find different 
criteria for the model parameters ruling the setting in of the polaronic regime. We also illustrate 
the common physical origin of polarons in the two models as arising from the simultaneous fulfilling 
of two conditions: A sizable ionic displacement and a lattice deformation energy gain larger that 
the loss in the bare electron kinetic energy. 



I. INTRODUCTION 



Polarons have been observed in various materials like 
Ti 4 C>7, Ni x V2C>5 and tungsten oxides WC>3_ X . More 
recently polaronic features have been detected in opti- 
cal spectroscopy in the liahtly doped insulating phase of 
cuprate superconductorsoa. 

Electrons acquire a polaronic character in the presence 
of a sufficiently strong electron-phonon {e-ph) coupling 
when they displace the ions around them and move car- 
rying along the lattice deformation. Being accompanied 
by the much heavier lattice degrees of freedom, in the po- 
laronic regime the carriers acquire large effective masses 
and, in some cases, may even be trapped in the poten- 
tial well arising from the ionic displacement that they 
created. 

Polarons are usually classified as being large or small 
polarons depending on whether the ionic deformation is 
spread over several sites or involves one single site. In 
this paper we confine ourselves to the investigation of 
models with short-range e-ph interaction, where, apart 
from a more or less-jrarrow crossover region only small 
polarons are formecB. _ 

In particular we investigate the Holstein modelQ and 
the Su-Schrieffer-Heeger (SSH) modefl. Due to their 
relative simplicity these models are definitely the most 
frequently considered models for electrons and phonons 
interacting via a short-range potential and may well be 
taken as a suitable paradigmatic basis for investigating 
the physics of strongly interacting e-ph systems. 

The concepts underlying polaron theory in these mod- 
els are long-standing and have found various theoreti- 
cal substantiations over the last decades. However, the 
strong-coupling nature of the polaronic state does not 
allow reliable analytic approaches in the intermediate 
crossover region, where one should sit in order to quanti- 
tatively investigate the conditions for polaron formation. 
On the other hand the multiphononic essence of polarons 
makes it difficult to approach the strong-coupling regime 
from the numeric point of view. A landmark in this con- 



text was provided by quantum Monte Carlo calculations 
in Ref. ^ where an interpolation formula was presented 
describing the critical e-ph coupling leading to polaron 
formation in the Holstein model. However, no detailed 
analysis and distinction between the adiabatic and anti- 
adiabatic regimes was carried out there and the interpo- 
lation formula was presented in a rather empirical way. 
More recently the discovery of polarons in the insulating 
phases of the high temperature superconductors has trig- 
gered numerical exact diagoaalization analyses on models 
with strong e-e interactionstjla. 

Despite the renewed interest, the literature still lacks 
a unifying picture providing a clear physical and formal 
understanding of the phenomenon of polaron formation 
both for the simple dilute case of one single polaron in an 
empty lattice and for the case of many interacting elec- 
trons. The main goal of this paper is to investigate the 
Holstein model and the SSH model in order to provide 
such a picture for the single particle case. 

In the Holstein model 
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the first term is proportional to the nearest-neighbor hop- 
ping integral t, which we will take as our unit of energy, 
giving rise to a d-dimensional tight-binding band struc- 
ture of the form E(k) — —2t J2t=i cos(k u ) (for simplicity 
we assume a cubic lattice with unit lattice spacing and 
h = c = 1) accounting for the kinetic energy of the free 
electrons. Since we will restrict ourselves to the single 
electron case we will not consider electron spin indices 
throughout this paper. 

A dispersionless (oj(q) — luq) Einstein phonon is cre- 
ated by the field a\ and is coupled to the local electronic 
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density in the Holstein model and to the covalent bond 
variable c\c i+ i + cJ +1 Cj in the SSH model. For the Hol- 
stein model, the coupling arises from the dependence of 
the local atomic energy (i.e. the Madelung energy) on the 
ionic position. This coupling is relevant when the screen- 
ing of the Madelung potential is poor, and is believed 
to be non-negligible in the superconducting cupratesll3. 
On the other hand the covalent e-ph coupling in the SSH 
model is due to the dependence of the hopping integral 
on the relative distance between two adjacent ionsEl. No- 
tice that our SSH model differs from the conventional 
one in having optical (instead of acoustic) phonons like 
the Holstein model: To clarify more easily the common 
mechanisms underlying the polaron formation in the two 
models wc avoided unnecessary differences between them, 
thus focusing on the role played by the different e-ph cou- 
plings. Moreover, despite the completely different origin 
of the e-ph coupling in the two models, we choose the 
same notation g to emphasize the generic character of 
the physical processes that we are going to present. 

Before addressing the problem of the single polaron for- 
mation in a more formal way within the above models, 
we first would like to provide simple and intuitive argu- 
ments. As mentioned above the setting in of a polaronic 
regime is characterized by both a lattice deformation en- 
ergy gain larger than the loss of bare kinetic energy and 
a sizable local displacement of the ionic positions, giving 
rise to a strong reduction of the effective hopping matrix 
element. 

For the Holstein model, these effects are directly re- 
lated to two parameters which are often introduced in 
this field: A = g 2 /(2dtLu ) and a = g/ujQ. The former 
has a twofold meaning. In fact it can represent the ef- 
fective phonon-mediated e-e coupling introduced, e.g. in 
Fermi liquid theory or in the traditional BCS theory of 
superconductivity. Alternatively and more interestingly 
in the present context, A may represent the ratio between 
the polaronic binding energy E p = —g 2 /ujQ (see below) 
in the strong coupling limit and the bare average kinetic 
energy of the electrons of the order of half the band- 
width (~ — 2td). Notice that the bare hopping t has to 
be used here. In fact this is of the order of the kinetic 
energy actually lost when the polaron is formed. Then 
the value of A determines the convenience for the sys- 
tem to give up the kinetic energy gain arising from the 
hopping to gain the lattice deformation energy induced 
by the local e-ph potential. On the other hand, as is 
clear from a standard Lang-Firsov transformation of the 
Holstein modeltll, a represents half the ionic displace- 
ment in units of (2Mu>q) where M is the ionic mass. 
The simultaneous occurrence of the two conditions A > 1 
and a > 1 is needed to characterize and to determine 
the polaron formation in the Holstein model. Thus one 
can immediately recognize from the definition of A and a 
that a crucial role is played by the adiabatic ratio u>o/t. 
If uja/t is small, the condition for a large A = a 2 cuo/2td 
is more difficult to realize than a > 1, and polaron for- 



mation will be determined by the more restrictive A > 1 
condition. The opposite is true when the system is in 
the anti-adiabatic regime u>o/t > 1. This intuitive argu- 
ment was already implicit in the interpolation formula 
Eq. (4.1) in Ref. || once this is expanded in the two op- 
posite limits LOo/t^> 1 and u>o/t -C 1. 

For the SSH model it will be shown that the value of A 
determines both the the tendency towards the localization 
of the electron and the suppression of the hopping inte- 
gral associated with the lattice distortion, regardless the 
value of the adiabatic ratio ui^/t. Therefore A will be the 
relevant parameter for the description of the system for 
any value of uo/t. This intuitive argument will be sub- 
stantiated below by analytic calculations and numerical 
exact diagonalization of small clusters. 



II. A FIRST INSIGHT FROM PERTURB ATIVE 
CALCULATIONS 

The above arguments may be made more formally pre- 
cise within a perturbative analytic calculation in the limit 
of small e-ph coupling (g <C t,LUo). In this case we eval- 
uated the second order correction to the electronic self- 
energy represented in the diagram of Fig. 1. 



FIG. 1: Lowest-order self-energy diagram of the electron 
propagator. The solid line is the bare electron Green func- 
tion; The dashed line is the bare phonon propagator. The 
dots represent the e-ph coupling. 

To explore both the adiabatic and the anti-adiabatic 
regimes, we cannot apply Migdal's theorem to discard 
vertex corrections: Our approximation only relies on the 
small value of g. This also allows for the simplification of 
using the bare electronic Green function instead of car- 
rying out a self-consistent evaluation including the full 
Green function inside the self-energy diagram. 

For the sake of simplicity and to allow for a more direct 
comparison with the numerical results on small clusters, 
we only deal with one-dimensional systems. Neverthe- 
less, we explicitly checked in the easier case of the Hol- 
stein model that our analytic results apply to multidi- 
mensional cases as well. 

For our purposes we only need the perturbative cor- 
rections to the effective mass 
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and to the ground state given by the solution of the equa- 
tion 



e - -ReS(fc = 0, w) = 0. 



(4) 



It is quite obvious that the polaronic regime cannot be 
attained within our lowest-order perturbative approach: 
The cloud dressing the electrons in a polaronic excitation 
involves multiphononic processes, which are not included 
in the diagram of Fig. 1. Nevertheless valuable indi- 
cations on the beginning of the polaronic crossover can 
be extracted from the above expressions. In particular, 
one can determine the parameters for which the effective 
mass starts to grow [(m* — m)/m ~ 1], also marking the 
region where perturbation theory is no longer applicable. 



A. The Holstein model 

For the Holstein model, in agreement with Ref. [l^ we 
obtained 
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where the real part of the square root has the same sign as 
(lo — wo)- Notice that the self-energy is momentum inde- 
pendent because the bare e-ph vertex is also momentum 
independent. This feature allows the introduction of the 
density of states in the momentum integrals, thus leading 
to a straightforward extension of our results above one 
dimension. 

In d = 1 we find 
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By expanding the perturbative correction in Eq.(j6|) i-u, 
the two opposite, adiabatic (ujq <§; t) and anti-adiabaticEJ 
(loq ^> t) limits, we get 
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In agreement with the intuitive arguments presented at 
the end of the previous section, the anti-adiabatic result 
shows that the mass enhancement is driven by the con- 
dition a > 1, which is more restrictive than A > 1. This 
finding also matches well with the result one would ob- 
tain from a standard Lang-Firsov transformation^ of the 
Holstein model. As is well known, this transformation 



changes the free electron basis into a polaron basis and 
shifts the ionic equilibrium position by a. By suitably 
choosing a — g/oJo the coupling between the phonons and 
the local electron density is eliminated and the fermion- 
phonon interaction is instead displaced into the kinetic 
energy term. In the anti-adiabatic limit, where t is small, 
this scheme is more convenient (the interaction is put in 
the small part of the Hamiltonian) and, averaging the 
Hamiltonian on the phononic vacuum one finds a bare 
polaron binding energy E p — — a 2 uiQ and an exponen- 
tial reduction of the kinetic energy t — > t* = £cxp(— a 2 ). 
This standard result can easily be connected to our per- 
turbative result (JsJ) by noticing that, for small g, a is also 
small so that t* ~ t(l — a 2 ). 

On the other hand, as expected, the adiabatic result 
for the mass enhancement is proportional to A. However, 
this contribution is also proportional to y/t/ujQ, which is 
large in this limit and this also seems to contrast the usual 
perturbative calculations within Fermi liquid theoryliil, 
which only predict a correction of order A. This feature 
arises from the singular density of states at the bottom 
of a one-dimensional band, as can be checked by con- 
sidering a finite density of electrons with a finite Fermi 
energy fj, away from the bottom of the band. In this lat- 
te r case the coefficient of the mass correction becomes 
y/t] L*p + A*) an d is no longer singular in the adiabatic 
limitlij. This specific, non-generic result is the price that 
we have to pay in order to take advantage of the sim- 
pler analytic treatment in one dimension, but it does not 
hide the important finding that only the A parameter 
rules the polaron formation when luq < t. Moreover it is 
worth noting that our simple, lowest-order perturbative 
calculation gives_a strong indication that one electron in 
a fully adiabaticE£l one-dimensional Holstein lattice is lo- 
calized (m* -r*,oo). In light of our calculation, this well 
known resultllZl can easily be attributed to the singular- 
density of statesE3. 



B. The SSH model 

The simplifying feature of a momentum independent 
bare e-ph vertex is no longer present in the SSH model, 
where the bare e-ph vertex in d — 1 associated with scat- 
tering of an electron from a k to a k + q Bloch state has 
the following form 



gk,k+ q = 2ig [sin(k + q) - sin(fc)] 



(10) 



Consequently, once a frequency integration is carried out, 
the self-energy is given by 



S(fe,w) = 4g 2 /N 
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x [sin 2 (A: + q) + sin 2 (fc) - 2sin(fc)sin(fc + q)] . (11) 

To obtain the desired physical quantities for a single elec- 
tron in the bottom of the band, we only need to evaluate 
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ReYj(k — 0,o;) and its derivatives in the k = state and 
for uj = E(k = 0) = -2t, finding 
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Also in this case we evaluate the mass correction in the 
adiabatic and in the anti-adiabatic limits 
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Contrary to the Holstein model case, the electrons in the 
SSH model in the fully adiabatic limit are completely free 
(to* = to). This is due to the vanishing of the bare e- 
ph vertex for small transferred momenta q [cf. Eq.QLB)], 
which overcompensates the divergent density of states in 
the integral of Eq. ( pi] ) . Physically this effect arises be- 
cause many electronic states lie close to the k — point, 
but the phonon-mediated scattering between them occurs 
at low momentum transfer, which is less effective in the 
SSH model, where phonons decouple from the electrons 
in the long- wavelength limit. 

An important difference between the Holstein and the 
SSH models is also present in the anti-adiabatic limit, 
where the mass enhancement in the SSH model [Eq.(15)] 
involves both a 2 and A. In this limit A ^> a 2 and there- 
fore the mass correction is dominated by A also in the 
anti-adiabatic regime. This difference arises from the 
different localization mechanism taking place in the SSH 
model with respect to the Holstein model. In light of this 
specific mechanism one can interpret the dependence of 
the mass enhancement on the parameter A in both the 
adiabatic and anti-adiabatic regimes. However, we pre- 
fer to postpone the discussion of this issue until the full 
numerical analysis of the SSH model is presented in the 
second part of the next section. This analysis will provide 
a clear substantiation of the above perturbative indica- 
tions and a natural interpretation will be given of the 
criterion for polaron formation. 



phonons per site (and check the convergence of the re- 
sults by varying the phonon number). Due to the huge 
enlargement of the Hilbert space induced by the pres- 
ence of the lattice degrees of freedom, we have only-been 
able to investigate small clusters up to four sitesea. In 
such small clusters finite size effects are obviously rele- 
vant. However, we checked that as far as the criterion 
for polaron formation is concerned, our results are rather 
insensitive of the boundary conditions and no qualitative 
changes occur in passing from three- to four-site lattices: 
In the short range models considered here, polaron for- 
mation is a local, high-energy phenomenon. 

The numerical calculations for the Holstein model have 
been performed using a slightly modified version of the 
Holstein Hamiltonian (|l|), in which the phonon displace- 
ment operator is coupled to the local electron density 
fluctuationaZi: 
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where rii = c\ci is the number operator for electrons 
on site i and (rii) is the mean value of the same quan- 
tity. This choice removes the trivial coupling between the 
mean electron density and the zero momentum phonon 
mode, and thus allows for better convergence as far as 
the number of phonons is concerned, without affecting 
physically relevant quantities. 
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III. EXACT DIAGONALIZATION ANALYSIS 

In order to found the above scheme for the single po- 
laron formation on more solid ground, we performed ex- 
act numerical calculations on small-clusters by means 
of the Lanczos algorithm. As usual! 2 !, we truncate the 
phononic Hilbert space so as to include only a finite 
number of phonons per lattice site. To reliably explore 
the strong-coupling regimes, we had to include up to 50 



FIG. 2: Local-density-displacement correlation function 
for the Holstein model and one electron in a four-site lat- 
tice with periodic boundary conditions. Lines are labelled by 

LUo/t. 

To extract information on the values of g at which 
the polaron crossover begins we analyzed the correlation 
function between the electronic density on a site i and 
the ionic displacement on the site i + 8 
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Xi,s = (4> \^2 c L c i<? ( a i+s + a l+s) IM (17) 

a 

where |</>o) is the ground state of the system. 

For the Holstein model we report in Fig. 2 the behav- 
ior of the local-density-displacement correlation function 
Xi,s=o f° r different values of the adiabatic parameter u)o/t 
as a function of the e-ph coupling constant. The calcu- 
lation was performed on a four-site lattice with periodic 
boundary conditions. As is well known, all physical quan- 
tities are smooth functions of the e-ph coupling. Never- 
theless it is evident that the adiabatic regime is character- 
ized by a rather sharp crossover, whereas the increase of 
Xi,<5=o is much slower when u>o/t > 1. By calculating non- 
local correlation functions with S ^ 0, we also checked 
that, in the Holstein model in one dimension, apart from 
the crossover region, polarons are always small. In fact 
the increase of the local Xi,0 is always accompanied in 
the strong-coupling regime by the decrease of the non- 
local correlation functions, showing that the polaron is 
so narrow that the presence of a fermion on a site is un- 
correlated with the ionic displacements on neighboring 
sites. 

Since polaron formation is a crossover without sym- 
metry changes between two phases, some arbitrariness is 
unavoidable in defining a criterion separating the free- 
electron and the polaronic regimes. In particular we 
choose the critical g from the point of maximum slope of 
the local-density-displacement correlation function Xi,o- 
We checked that different criteria (like, e.g., the maxi- 
mum of the nearest-neighbor-density-displacement cor- 
relation function Xi,i) provide the same qualitative re- 
sults. 
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FIG. 3: Phase diagram for the Holstein-model with one 
electron in a four-site lattice with periodic boundary condi- 
tions: The crossover region is shaded, while the solid line is 
the critical value defined in the text as a function of A (a) and 
of a (b). 

From the critical values of g we calculated the corre- 
sponding values of A obtaining the phase diagram of Fig. 
3(a). The critical A as a function of wo/t is indicated with 
a solid line. We also show the crossover region, defined 
as the range of parameters for which Xi,o significantly 
changes (we calculated the second derivative of Xi,o to 
numerically evidentiate this region). 

This clearly shows that at small phonon frequencies 
the criterion for having a polaronic regime is A = A c ~ 1 , 
whereas at larger phonon frequencies one obtains uio/t ~ 
2A C , which implies a c w 1. In Fig. 3(b) the same in- 
formation of Fig. 3(a) is translated in terms of a. In 
this way it is made clear that a c is constant in the anti- 
adiabatic regime uj > t. Therefore, the exact numerical 
solution on a small cluster confirms the intuitive argu- 
ments as well as the perturbative calculation!^. 

The same analysis has been carried out for the SSH 
model. While the general scheme is the same, some addi- 
tional care has to be used here to check that the adopted 
values of g do not invert the sign of the hopping parame- 
ter or of the kinetic energy. In the SSH model, an electron 
tends to shrink a bond increasing the effective hopping 
between two sites. At the same time the neighboring 
bonds are stretched and the hopping between the two oc- 
cupied sites and the surrounding ones is reduced resulting 
in a tendency towards localization. Due to translational 
invariance this hopping reduction is translated into an 
effective reduction of the quasiparticle bandwidth (i.e. 
an enhancement of the effective mass). Eventually the 
hopping between the two sites and the rest of the lattice 
vanishes and may even change sign. This pathological sit- 
uation is a well known feature of the SSH model, which in 
real systems never occurs due to higher order corrections 
to the expansion of the hopping parameter t in terms of 
the ionic displacement. For all couplings where we find 
polaron formation in the SSH model, we checked that 
these pathologies do not occur. We also notice that, by 
increasing the phonon frequency, the effective hopping 
is relatively less affected, so that larger values of A can 
be reached before the zero-hopping pathology is found. 
As a consequence, the region with substantial polaronic 
character is enlarged. 

The results for the local density-ionic displacement cor- 
relation functions do not qualitatively differ from those 
of the Holstein model. Again one finds that for finite 
phonon frequencies physical quantities are smooth func- 
tions of the e-ph coupling, and that the polaron forma- 
tion is a rather fast crossover when cuq < t, whereas the 
crossover is slower for larger loq. In the present case the 
phase diagram is reported in Fig. 4. 

This diagram is evaluated making use of the nearest- 
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neighbour correlation function, which is an increasing 
function of g in weak-coupling and it is decreasing in 
strong-coupling. The maximum of this function will then 
be the critical value for polaron formation. The crossover 
region is estimated by the "width" of the same function, 
whereas the pathological region of parameters is associ- 
ated with a negative value of the same quantity, as a 
consequence of the unphysical negative value of the ef- 
fective hopping matrix element. 

As expected from the perturbative calculation, we find 
that the polaronic regime is determined by the condition 
A c =const. both in the adiabatic and the anti-adiabatic 
regimesEil. The specific mechanism of hopping reduction 
giving rise to localization in the SSH model, accounts for 
this difference with respect to the Holstein model. 
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FIG. 4: Phase diagram for the SSH-model with one elec- 
tron in a four-site lattice with open boundary conditions: The 
solid line is the critical value of A for the polaronic crossover; 
The crossover region is shaded and its boundary is indicated 
by the dashed lines. The dot-dashed line is the boundary of 
the pathological region. 

For clarity we begia this discussion by considering the 
fully adiabatic limitll3. In this case, at first order the 
hopping is reduced by the stretching of the bonds by 
t* = t — gu, where u is the (dimensionless) value of the 
bond length variation in units of (2Mwo)~ 5 (notice that 
for luq — > 0, g has to vanish but gu stays finite)ll3. Since 
it is possible to show that also in the SSH model the lat- 
tice displacement [in units of (2Mwo)" 5 ] is proportional 
to g/u>o, one finds 



r = t-gu = t 1-7 
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In the strong coupling regime the electron is localized 
on a single bond and the many-site model becomes equiv- 
alent to a two-site cluster. Then, for such a small sys- 



tem, an analytic solution is possible giving 7 = 4 for 
the constant appearing in Eq. (|l8|) . Thus A determines 
the reduction of the effective hopping when the lattice 
displacement is sizable. 

As for the Holstein model case, another condition 
needs to be satisfied in order to have a polaronic regime: 
The energy gain due to the polaron formation, i.e. the 
typical energy of a self-trapped carrier has to be larger 
than the loss in bare kinetic energy associated with the 
self trapping. Since in the strong coupling regime the 
electron localizes on a single bond we can limit ourselves 
to a two site-cluster in order to evaluate the polaronic 
binding energy. The solution of the two-site cluster shows 
that the polaron energy for the SSH model still contains a 
contribution from the free electron hopping, arising from 
the derealization of the electron between the two sites 
of the bond. Then the ground state energy is given by 
Eq = —t — 2g 2 /ujq. We divide this energy by the free 
electron energy — 2t to obtain the ratio of the energy gain 
associated with polaron formation to the energy-loss as- 
sociated with the decrease of electronic mobilityc3. If we 
explicitly evaluate the range of A values for which 



(-t - 2g 2 /LU )/(-2t) > 1, 



(19) 



we readily obtain A > 0.25. This value coincides with 
the value at which the hopping matrix element vanishes 
according to Eq.(|I§|) and to the adiabatic limit of the 
parameter 7 = 4. This implies that when luq — ► 0, 
the system will have no energetic advantage in localizing 
the electron on a bond, unless the pathological condition 
t* = is reached. According to the physical idea that 
both a sizable lattice displacement and an energy gain 
from deformation larger than the kinetic energy loss are 
required to realize a polaronic state, one should not ex- 
pect polarons in the adiabatic limit of the SSH model. 
Indeed, we carried out the exact diagonalization of large 
clusters (100 sites) in the extreme adiabatic limit finding 
that the SSH model does not present any marked pola- 
ronic behavior for couplings smaller than the "patholog- 
ical" g's at which the hopping changes sign. 

For finite phonon frequencies this picture is modified 
by the lattice dynamics. The numerical study shows that 
the ground state energy is not strongly effected by the 
lattice dynamics: regardless of the value of wrj/f, A larger 
than 0.25 remains the condition to obtain an energetic 
advantage from localization. On the other hand the ef- 
fective hopping matrix element is less severely reduced 
by the coupling to the lattice fluctuations and the value 
of A for which the effective hopping becomes zero in- 
creases with u>o/t (from Fig. 4, one sees that 7 w 2 for 
LUo/t Ri 20) Therefore, for finite luq, it is always possible 
to find a regime where the lattice deformation becomes 
energetically favorable and a substantial lattice displace- 
ment (i.e. hopping reduction) is present without having 
a non-physical vanishing of the hopping. In this region, 
which is larger for large phonon frequencies (see Fig. 4), 
the electron has a polaronic character for values of A 
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larger than A c w 0.2. Finite size effects easily account 
for the small quantitative discrepancy between this value 
and fully adiabatic estimate A c w 0.25. 

The above result shows that A determines both the re- 
duction of the hopping integral associated with the lattice 
distortion and the tendency towards localization driven 
by the energetic advantage in deforming the lattice. It is 
then natural to consider A as the relevant parameter for 
polaron formation regardless of the value of the adiabatic 
ratio ojo/t. Notice that this finding was also suggested by 
the perturbative result showing that the main cor- 
rections to the effective mass are proportional to A both 
in the adiabatic and in the anti-adiabatic regime. 



models. 
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IV. CONCLUSIONS 

In the present paper we addressed the issue of polaron 
formation in lattice models with extreme short-range e- 
ph interactions, the Holstein and the SSH models. Our 
work was devoted to settle in a definitive way the case of 
a single polaron in a dynamical lattice. Despite compre- 
hensive and repeated investigations, this topic still lacked 
a clear systematic conclusion, leading to sometimes con- 
tradictory and confusing statements being reported in 
the literature. In particular we clarified, both from an- 
alytic qualitative arguments and from numerical exact 
calculations, that A or a are not by themselves indepen- 
dent parameters which determine the free-electron or the 
polaronic regimes in the e-ph models. Indeed we showed 
for the Holstein model that both conditions, A > A c w 1 
and a > a c w 1 have to be satisfied in order to realize 
both the kinetic energy reduction and the sizable lattice 
displacement which characterize the polaronic state. De- 
pending on the adiabatic ratio loq /t, the condition for the 
polaronic regime is determined by A > A c when loq < t 
and by a > a c when lvq < t. 

Comparing our findings with the results of a dynamical 
mean field theory calculation, which is exact in the limit 
of infinite connectivity!^!, we find substantial agreement 
as far as the value of the parameters ruling the single- 
polaron formation in the Holstein model is concerned. 
This clearly indicates that the same physical picture ex- 
tracted here from the numerical calculation in small (one- 
dimensional) clusters holds for infinite systems in higher 
dimensions as well. This "universal" behavior is a natural 
consequence of the local character of the small polarons 
in the Holstein model. 

On the other hand, perturbative and numerical calcu- 
lations for the SSH model lead to the condition A > A c ir- 
respective of the adiabatic ratio. This result, apparently 
contrasting with the corresponding criterion in the Hol- 
stein model, may be very naturally understood in terms 
of a physical argument based on the kinetic energy re- 
duction and the ionic deformation. Therefore, although 
the final criteria are different, the same physical picture 
underlies the formation of a single polaron in the two 
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